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reviewed  extensively[l-4].The  magnetic  susceptibility  of  dimeric  complex- 
es has  also  been  the  subject  of  several  reviews[3,5-7],  and  a number  of 
n.m.r.  contact  shift  studies[8]  of  dimeric  complexes  have  been  reported. 
The  Interpretation  of  all  such  studies  rests  upon  the  determination  of 
the  parameters  contained  in  the  spin  Hamiltonian  of  the  pair.  The  eval- 
uation of  matrix  elements  then  is  a primary  concern,  and  there  are  basi- 
cally two  ways  In  which  to  approach  the  problem.  In  the  first  method 
one  takes  advantage  of  the  fact  that  the  spins  are  coupled  by  taking 
linear  combinations  of  single  ion  state  vectors,  so  that  the 

resulting  wave  functions  are  eigenfunctions  of  the  total  spin  operators 
o 

z and  1(1^  = but  treats  the  problem  completely  in  terms  of 

single  spin  operators.  In  the  second  method  one  starts  with  the  coupled 
representation  state  vectors,  |SM>,  and  transforms  the  Hamiltonian  so 
that  it  contains  only  total  spin  operators.  The  generality  of  this 


method  provides  insight  that  is  lacking  in  the  first  approach  to  the  pair 
problem  and  also  reduces  the  number  of  operator  expressions  required 
to  evaluate  the  Hamiltonian  matrix.  To  our  knowledge,  a full  treatment 
of  the  general  case  using  the  "total-spin"  Hamiltonian  has  not  been 
reported.  Thus,  it  is  the  purpose  of  this  paper,  to  report  the  results 
obtained  by  the  latter  method  through  the  application  of  the  Wigner- 


Eckhart  theorem. 
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I.  The  Pair  Hamiltonian. 


We  assume  that  the  pair  Hamiltonian  can  be  expressed  as  the  sum  of 
the  two  single  Ion  spin  Hamiltonian  plus  terms  representing  Interactions 
between  the  two  Ions  of  the  pair.  This  implies  that  the  interaction 
terms  are  small  compared  to  the  crystal  field  splittings  and  will  thus 
have  a negligible  effect  on  the  single  ion  tensors.  The  general  form 
of  the  pair  Hamiltonian  can  then  be  written  as  follows: 

H — ff-j+^2  ~ 2JS^ '$2  + $1  * D g*-d^2*  0) 


where 


"i  * «-§rSt  ♦ vsrsi  * VMi 


Eq.  (1)  assumes  nothing  about  the  relative  orientation  or  form  of  the 
various  tensors,  and  can  contain  anisotropic  exchange  terms  as 
well  as  the  dipole-dipole  contribution.  The  usual  procedure  for  eval- 
uating the  matrix  of  (1)  Is  to  expand  the  total  spin  state  vectors 
|S/i>,  in  terms  of  the  single  spin  state  vectors,  as  follows: 


| SM>  = l C(S1S2S;M1M2M)|S1S2M1M2>,  (2) 

m1m2 

where  is  the  Clebsh-Gordan  coefficient,  for  which  we  have 

adopted  Rose's [9]  notation.  Once  the  wave  functions  are  in  this  form  one 
can  use  (1)  directly  to  evaluate  matrix  elements.  While  this  process 
is  straightforward,  it  is  quite  tedious,  and  the  expansion  in  (2)  must 
be  found  for  each  pair  system  characterized  by  different  S1  and  S2. 

It  Is  possible  to  avoid  the  expansion  in  (2)  by  evaluating  the  matrix 
of  (1)  In  the  coupled  representation,  and  to  facilitate  this  process 
we  have  rewritten  (1)  as  follows: 
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H * - S1(S1+1)-S2(S2+1)} 

+ B/2{H-^-Z  + H-§5  -A}+  1/2{Z*A1*I1  + A-A^} 

+ 1/2{E-A2-I2-A-A2'I2}+  V4{l--Defd'£-A.-Pefd-I} 

+ l/4{r-D  -Z+A-D  ’A}  + 1/2Z'D/A,  (3) 

where  Z * S-j+Sg,  - = -l“-2’  -o  = -l+-2*  ^5 = -l“-2*  etc'  state 

vectors  |S,M>  are  eigenfunctions  of  the  familiar  total  spin  operators 

z*z  and  Z2,  and  matrix  elements  of  the  form  <S' M1  | z . | SM> , (1  = x,y,z), 
vanish  for  S'/S.  Terms  containing  components  of  the  less  familiar  vec- 
tor operator  a,  however,  follow  the  selection  rule  S'  = S,  S+l,  S+2. 

The  matrix  elements  of  terms  containing  (i  = x.y.z)  can  also  be 
evaluated  In  the  coupled  representation.  Much  of  what  will  follow  then, 
will  deal  with  finding  the  matrix  elements  of  the  a^  and  in  particular 
we  shall  attempt  to  show  two  important  results: 

1.  For  the  special  case  of  matrix  elements  within  a given  S mani- 
fold (<S'M'|tf|SM  >,  where  S = S'),  any  operator  expression  containing 
A.,  can  be  replaced  by  the  corresponding  expression  containing  Z^  if  a 
proportionality  constant  is  applied.  This  result  is  particularly  impor- 
tant since  if  -2JS^*S2  is  much  larger  than  the  other  terms  in  B ( a case 
often  encountered  experimentally),  these  are  the  only  matrix  elements 
one  needs  to  evaluate. 

2.  For  matrix  elements  between  spin  manifolds  (S7S),  the  matrix 
elements  of  a^  can  be  evaluated  by  general  formulas  that  are  simple  to 
use,  and  make  possible  the  identification  of  zero  matrix  elements  with- 
out tedious  algebraic  manipulations. 
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II . Irreducible  Tensor  Operators  and  The  Wigner-Eckhart  Theorem. 

For  much  of  what  will  follow.  It  will  be  convlenlent  to  Introduce  j 

the  Irreducible  tensor  operators.  These  operators  are  defined  by 

j 

Racah  through  the  following  commutation  relationships: 


CV  Tlm(x)]  ‘ wTZm(x)-  (4) 

[T+J^x)]  = {i(  z+l  )-m(nH-l ) >1/2Tz+1  #m(x) 
where  T^m(x)  Is  the  m^-  component  of  an  Irreducible  tensor  operator  of 
rank  l and  operator  variable  x.  The  operator  variables  of  interest 
here  are  x = e,a,S-|,  and  Sg.  It  is  also  possible  to  construct  irreduc- 
ible tensor  operators  of  mixed  operator  variable  and  for  the  present 
case  we  shall  need  two  such  operators  of  this  type,  T2*/SlV  and  T2m 
(ea).  In  Table  I we  have  collected  all  of  the  irreducible  tensor  opera- 
tors we  shall  have  need  of,  and  It  Is  easy  to  show  that  they  do  Indeed 
satisfy  the  commutation  relationships  In  Eq.  (4).  It  can  also  be  seen 
from  Table  I that  the  Irreducible  tensor  operators  are  nothing  more  than 
normalized  linear  combinations  of  ordinary  cartesian  vector  operators. 

A theorem  involving  irreducible  tensor  operators  that  we  shall  use 
extensively  is  the  Wigner- Eckert  Theorem,  which  can  be  expressed  as 
follows: 

<S'M’  lTZm<x)lSM>  = C(SJS’ ;MnM' )<S' | |Tz(x)||S>.  (5) 

Thus,  the  Wlgner-Eckart  theorem  states  that  the  matrix  elements  of 
an  irreducible  tensor  operator  can  be  decomposed  Into  the  product  of 
a Clebsch-Gordan  coefficient  and  a quantity  commonly  referred  to  as  a 
reduced  matrix  element.  The  Clebsch-Gordan  coefficient  contains  all  of 
the  directional  or  m dependence  and  is  conpletely  independent  of  the 
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operator  variable,  x.  The  reduced  matrix  element  contains  purely  dyna- 
mlcal  Information  and  Is,  as  the  notation  Implies,, completely  Independent 
of  If  , M and  m.  It  Is  obvious  that  (5)  can  be  extended  to  Include  sums 
of  Irreducible  tensor  operators  of  the  same  rank  so  that  we  may  write 
formally, 

<S'M' |ETlm(x) |SM>*<S* | |Tj(x) 1 1 S>EC(S  ZS ' ;MmM' ) . (6) 


m m 

The  Importance  of  Eq.  (6)  lies  in  the  fact  that  all  of  the  operator 


expressions  that  arise  from  the  pair  Hamiltonian  (Eq.  (3))  can  be 
expressed  as  linear  combinations  of  Irreducible  tensor  operators.  To 
show  this  Is  a relatively  simple  matter  which  we  will  illustrate  for 


one  of  the  operators  In  (3),  a'Do-a. 

0 0 


a-D^-a  = a-  C 

Li 


0^OXX  D 


° n 
0 ^ zz 


0 A 
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where  Dq  Is  in  its  diagonal  representation.  If  we  choose  some  non- 
diagonal representation  for  D . the  expression  for  a-D  'a  will  be  more 

—a  — — o — 

complicated,  but  it  will  still  prove  possible  to  express  It  as  a linear 
combination  of  irreducible  tensor  operators.  All  of  the  operators  In 
Eq.  (3)  can  be  treated  In  a similar  manner,  and  once  this  is  done  the 
Wigner-Eckart  theorm  (Eqs.  5 and  6)  can  be  used  to  take  matrix  elements. 
We  have  done  this  for  all  operators  containing  a,  and  collected  the 
results  in  Table  II.  Many  of  the  operator  expressions  in  Table  II 
will  occur  only  if  the  interaction  tensors  are  non-diagonal.  In  order 
to  obtain  the  equivalent  expression  for  operators  involving  z one 
simply  replaces  A by  z everywhere  in  the  expression. 

Table  II  shows  that  we  have  reduced  the  task  of  finding  the  matrix 
of  (3)  to  finding  sums  of  Clebsh-Gordan  coefficients  for  which  standard 
tables  are  available[10],  and  evaluating  reduced  matrix  elements.  The 

reduced  matrix  elements  will  be  considered  in  Section  III. 

An  additional  expression  of  interest  can  be  derived  by  applying 
Eq.  (6)  twice,  namely: 

<S'Vt  |tTZm(x)|SH>*{<S||Tz(x)||S>/<S'  ||Ti(z)||S>> 
m 

x<S'H'  | ZT^zJISM^ 
m 

However,  the  above  becomes  Indeterminate  for  S f S'  since  both 

<S'MlzTZm(z) |SM>  and  <S' | |Tz(z) | |S>  are  zero  in  this  case  and  thus 
m 

we  are  left  with  the  less  general  expression: 


<SM'|lTta(x)|SM> 

m 


PsIIVkHIs.-],  <sm' Iit,  (e)Ism>.  (7) 

pnrfernsd  ■ 


Eq.  (7)  states  that  within  a given  S manifold,  the  matrix  elements  of 
any  Irreducible  tensor  operator  of  operator  variable  x are  proportional 
to  those  of  the  corresponding  one  with  operator  variable  z.  The  pro- 
portionality constant  Is  the  ratio  of  reduced  matrix  elements  and  this 
will  be  different  in  general  for  different  S manifolds.  We  will  be 
interested  in  two  cases,  namely  x = a and  x = za,  and  (7)  implies  that 
for  S'=S  elements  we  shall  have  no  need  for  Table  II  since  we  can  make 
the  following  operator  substitutions: 

H*G6-A  = R^H'V-  ’ -^6  - = R' 

E-P6-A  = R"^S^Z-D6'Z. 

where  the  R^  are  reduced  matrix  element  ratios.  The  advantage  of 
this  is  that  the  matrix  elements  of  operators  containing  only  z are 
easily  evaluated  in  the  coupled  representation,  and  thus  once  the 
appropriate  reduced  matrix  element  ratios  are  found  we  can  evaluate 
any  S'=S  element  with  the  use  of  ordinary  operator  algebra. 

Ill . Reduced  Matrix  Elements 

In  order  to  evaluate  the  reduced  matrix  elements  in  Table  II  we 
again  use  Eq.  (5),  the  Wigner-Eckart  Theorm.  Then,  for  x = a we  have: 

<S'M' |T20(a)|SM>  * C(S2S',M0M^<S' | |T2(A)| |S> 

It  is  easily  shown,  however,  that 

<s'M'|t20(a)|sm>  = <s,mi|t20(s1)+t20(s2)-2t20(sts2)|sh> 

* <S*M1|2T20(S])+2T20(S2)-T20(Z)|SM>. 
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Applying  (5)  to  this  expression  yields 

<S’M,|T20(A)|SM>  = C(S2S,;MOH,){2<S'||T2(S1)||S>+2<S,||T2(S2)||S> 

- <S'||T2(e)||S>}. 

Comparing  the  two  results  we  find  that  we  have  an  expression  for  the 
reduced  matrix  elements  of  T2(a)  In  terms  of  those  for  T2(S^),  T2(S2), 
and  T2(e).  This  procedure  can  be  applied  to  the  other  reduced  matrix 
elements  In  Table  II  and  the  results  are  given  In  Eqs.  (8). 

<S * | |T1 (a) | |S>  = <S,||T1(S1)||S>  - <S'||T1(S2)||S>  (8a) 

<S'||T2(ea)||S>  * <S'||T2(S1)||S>  - <S'||T2(S2)||S>  (8b) 

<s* | |T2(E)| |S>  - 2<S* I |T2(S1 ) I |S>+2<S' I |T2(S2) | |S>-<S' I |T2(E)| JS>^ 

The  advantage  of  using  Eqs.  (8)  Instead  of  deriving  theoretical  expres- 
sions for  Tj(a)  and  T2(ea)  directly.  Is  that  general  expressions  for 
<S'I|TZ(S,)||  S>  have  already  been  derived  (see  for  Instance  Reference 
2.,  Rose).  For  example.  It  can  be  shown  that, 

<S'|JTZ(S1)||S>  * (-l)S2+Z'Sl'S  [(2S1+1)(2S+1)]1/2W(SS'S1S1;- 

JS2)<S1||Tz(S1)||S1>,  (9a) 

where  W(SS'S^S^;  ZS2)  Is  the  Racah  coefficient;  standard  tables  of 
Racah  coefficients  similar  to  those  for  Clebsh-Gordan  coefficients  are 
available[ll],  Eq.  (9a)  then  relates  the  matrix  of  T^(S1 ) in  the  coupled 
representation  to  Its  matrix  In  the  uncoupled  representation.  The 
expression  for  TZ(S2)  is: 


WIN 
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<S*  I |TZ(S2)|  |S>  ■ (-l)Sl+t'VS[(2S2+l)(2S+1)]1/2W(SS,S2S2;ZS1) 

<S2|  |Tz(S2)  I |S2>,  (9b) 

where  the  change  in  phase  from  (9a)  to  (9b)  can  be  thought  of  as  arising 
from  permuting  the  state  vectors  | S2S-j S>  to  their  usual  form  |S^S2S>. 

The  remaining  reduced  matrix  elements  that  must  be  considered  are 
<S*  | |Tj(z) | |S>,<S..  | |T^(S^ ) | |S •>  (i  = 1,2).  This  can  be  done  directly 
through  the  Wigner- Eckart  Theorm  as  we  shall  illustrate  for  Tjm(l). 

From  Eq.  (5)  we  have, 

<S'| |T] (E) | |S>  ^S'M1 |Tlm(E)|SM>/C(S2Sl;MmM'). 

Since  the  l.h.s.  must  be  independent  of  M',M  and  m,  we  choose  M'  = M, 
m = 0 and  obtain: 

<S'|T](Z)||S>  = <S'M|T1£)(z)|SM>  = <S'M|ZZ|SM> 

C ( S 1 S * ;M0M)  C(S1S' ;M0M) 

where  use  has  been  made  of  Table  I in  the  last  step.  Evaluating  the 
matrix  element  and  using  a table  of  Clebsh-Gordan  coefficients  we  have: 


<S* | | T] (e ) | |S> 


= VsM 


ff/[S(S+V)TT/2 


=6s,s[S(S+l)]1/2. 


The  same  procedure  can  be  used  to  evaluate  <S' | |T2(z)|  |S>,  <-S^  j |Tj(S-| ) | |S^>, 
etc. 

We  now  have  expressions  for  all  of  the  reduced  matrix  elements  of 
T^(S-j ) , Tj(S2)  and  T^i),  and  these  yield  explicit  algebraic  formulas 
through  the  use  of  standard  tables  of  Racah  and  Clabsh-Gordon  coeffic- 
ients; these  formulas  are  collected  in  Table  III.  Thus,  with  the  aid 
of  Eq.  (8)  and  Table  III  we  will  be  able  to  tabulate  explicit  express- 
ions for  all  of  the  matrix  elements  in  Table  II  and  also  to  evaluate 

f 

any  reduced  matrix  element  ratio  we  wish. 


1 
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IV.  The  Matrix  of  the  Pair  Hamiltonian 

From  the  foregoing  it  follows  that  a convenient  way  to  evaluate 
the  matrix  of  (3)  is  to  consider  it  as  two  separate  Hamiltonians,  one 
which  connects  S'  = S elements  and  one  for  S'^  S elements.  It  then 
follows  from  (7)  that  the  former  will  contain  only  Z as  an  electronic 
spin  operator  and  thus,  these  matrix  elements  can  be  evaluated  by  the 
usual  operator  methods  as  soon  as  the  proportional ity  constants  In  (7) 
are  known.  For  S’^S  elements  the  Spin  Hamiltonian  will  have  only  those 
spin  operators  containing  £ and  we  will  need  a separate  expression  for 
each  such  matrix  element.  The  remainder  of  this  paper  then  will  be 
aimed  at  illustrating  how  one  goes  about  using  the  formalism  of  the 
previous  sections,  and  we  shall  try  to  do  this  in  such  a way,  that  one 
may  use  the  formulas  without  becoming  enmeshed  In  the  details  that  led 
to  them. 

A.  S'=S  Elements 

We  begin  by  rewriting  (3)  entirely  in  terms  of  the  electronic  total 
spin  operator  £ according  to  (7)  as  explained  earlier: 


V=s  * -J<L*I-si(si+i)-s2(s2+D> 


(ID 


+(e/2)H-{6  +CG.1T 

— — o “o  — 


♦ (l/2)£‘(<'-C*)!Wd  * Me*.?**’1  • 

where  C and  C+  are  functions  of  S,  S1  and  Sg,  and  i,  Gg.  Gfi,  etc.  have 
been  defined  below  (3).  From  this  it  follows  that  for  each  of  the  spin 
multlplets  C and  C+  will  differ,  and  thus  one  is  using  a somewhat 
different  Hamiltonian  for  each  manifold  of  spin  S.  The  definitions  of 
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these  constants  In  terms  of  reduced  matrix  element  ratios  and  explicit 
algebraic  formulas  for  them  are  given  In  Table  IV.  Otherwise,  one  uses 
(11)  In  the  usual  manner,  with  coupled  representation  state  vectors, 

|SM>,  as  a basis  set. 

If  one  considers  the  case  where  J is  much  greater  than  the  other 
terms  In  (3)  we  can  consider  the  properties  of  the  pair  as  arising 
from  Isolated  spin  multiplets,  and  use  (11)  exclusively.  Experimentally, 
this  results  in  a rather  interesting  e.p.r.  spectrum;  each  multlplet 
can  be  thought  of  as  arising  from  an  effective  spin  Hamiltonian  of  the 
form: 

H = eH-G'Z  + E-A^Il  + E*A2*I2  + E-0»Z 
and  in  general,  these  "effective"  G,  A and  £ tensors  will  be  different 
for  each  spin  state.  These  tensors  are  not,  however,  independent  but 
are  related  to  the  f undemental  tensors  of  the  system  through  the  con- 
stants in  Eq.  (11).  For  S^Sg,  all  effective  £ and  A^  tensors  will  be 
equal  to  the  actual  G and  A tensors  for  all  spin  states  since  Table  IV 

reveals  that  C * 0 in  this  case.  An  interesting  example  of  a case  where 

all  spin  states  were  observed  by  e.p.r.  has  been  reported  for  Cr  pairs 

doped  in  the  spinel  lattice  MgAl^.  These  authors  were  aile  to  relate 

the  spectrum  of  each  spin  state  to  the  same  set  of  fundamental  tensors. 
Also,  they  were  able  to  observe  the  effect  of  the  tensor  unobscured 
by  the  other  possible  contributors  (D^  and  Dfi)  to  the  "effective"  D 
tensor.  Examination  of  Table  IV  and  Eq.  (11)  shows  that  for  S^=S2*3/2 
this  will  always  be  possible  since  for  the  S*2  state  spectrum  C+  = C_  * 0 
and  thus  will  be  the  only  fine  structure  term. 
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It  should  be  pointed  out  that  other  Hamiltonians  similar  to  (11) 
have  been  suggested  previously.  Owen[l]  has  given  a pair  Hamiltonian 

appropriate  for  the  special  case  - S2  and  all  tensors  sharing  the 
some  major  axes.  Choa[13]  has  derived  a more  general  Hamiltonian  similar 
to  (11),  which  in  our  notation  becomes: 

Hs  = -J{£*E-S1(S1+1)  - S2(S2+1)}+B/2H*{K1G1  + K2G2}-E 

+1/2E -{K1A1  •I1  + K2A2«I2}+r{l/2De+<J-l/2C]DeKj+C1D1  (12) 


By  definition,  Choat  constants  should  be  related  to  ours  as  follows: 

C = lyi  = 1 - K2,  C+  = l^U^+Cg).  (13) 

Using  (13)  and  the  expressions  below  (3)  (i.e.  6^=  G^  + G2,  §6  = §1  - G2, 
etc.)  one  can  show  that  (11)  and  (12)  are  In  fact  identidal.  Unfortunat- 
ely, comparing  Table  IV  with  the  algebraic  formulas  in  reference  13 
reveals  that  while  the  definitions  involving  K^,K2  and  C in  Eq.  (13) 
hold  those  involving  ,C2  and  C+  do  not.  The  error  is  apparently  In 
Choas  expressions  for  C-j  and  C2  since  with  these,  Eq.  (12)  gives  incorrect 
matrix  elements  as  compared  to  either  (1)  with  the  expansion  in  (2) 
as  a basis  set  or  (11)  with  the  |SM>  as  a basis.  The  latter  two  have 
given  identical  matrix  elements,  as  they  should,  for  all  cases  we  have 
checked.  Since, 


<S||T2(S,)||S> 
C1  * <5||T2(e)1|S->' 


<S||T2(S,)||S> 

and  C2 s <s||T2(E)||sr-  • 


the  correct  formula  for  these  constants  can  be  found  in  Table  III.  It 
is  also  worthwhile  mentioning  that  for  the  special  case  S-|=S2,  C+  1S 
identical  to  the  Owen-Judd  constant[l]  as  it  should  be. 


8.  S'^S  Elements 

As  discussed  earlier,  for  these  matrix  elements,  one  need  only 
consider  operators  involving  a since  those  not  containing  this  are  zero. 
Eq.  (3)  then  reduces  to: 


= &/2!T!v-  + - a-a2-i2} 

(14) 

+ l/4{A-Po*A+2r-p6-A-A‘De+<j-A}. 

The  matrix  elements  of  all  possible  operator  expressions  that  can  arise 

from  (14)  are  given  in  Table  V;  these  formulas  result  from  direct  appli- 
cation of  the  Wigner-Eckart  Theorem,  Eq.  (5).  Table  V then,  is  just  an 
explicit  form  of  Table  II,  and  thus,  the  expression  for  each  matrix 
element  is  the  product  of  two  quantities:  a Clebsch-Gordan  coefficient 
sum  and  a reduced  matrix  element.  In  using  (14)  one  should  keep  In  mind 
that  it  is  easy  to  identify  the  non- zero  matrix  elements  even  before 
Table  V is  consulted.  This  follows  from  two  considerations: 

1.  The  two  states  to  be  connected  must  satisfy  the  triangular 
relationship  with  the  rank  (the  l of  the  TZm(x))  of  the  connecting  opera- 
tor. In  (14)  we  have  only  first  rank  (H*Gfi-A,  a*A^*I^)  and  second  rank 
(a-Do-a,  S-Dg-A,  A*pe^j*A)  operators.  Thus,  the  former  can  connect  two 
states  only  if  S'  * S+  1 , and  the  latter  only  If  S'  = S+  1,  S + 2,  (dis- 
regarding of  course  S'  * S elements  which  have  already  been  considered). 
Furthermore,  it  turns  out,  that  z*p6 -a  connects  only  S'  ■ S + 1 elements 
even  though  it  can  be  expressed  as  sum  of  second  rank  Irreducible  tensor 
operators. 


■’* 
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2.  The  operator  expressions  Involving  A^(1  = x,  y,  z)  havz  pre- 
cisely the  same  "M"  dependence  as  their  Ei  counter  parts.  To  choose  a 

concrete  example,  3^2-a*a  and  3eza2-e* A connect  only  those  states  where 

2 2 2 

H'*M  just  like  their  familiar  3e z-e*e  analog.  Similarly  Ax-Ay  and 
E a - E a can  connect  only  those  states  where  M'  = M+2  as  Is  true  for 

xx  y y 

e2-e2.  Knowing  the  above  properties  of  (14)  greatly  reduces  the  number 
x y 

of  matrix  elements  one  need  consider  and  It  should  be  noted  that  most  of 
what  has  been  said  In  this  regard  follows  Immediately  from  the  Wlgner- 
Eckart  theorem. 

Examination  of  Table  V shows  that  we  have  given  all  of  the  explicit 
formulas  necessary  to  determine  the  unique  half  of  the  Hamiltonian  matrix 
for  S'/S  elements;  the  other  half  can  be  found  through  the  relationship, 
<S'M' |SM>e<SM|ff|S'M'>**  To  our  knowledge  this  Is  the  first  time  such 
formulas  have  ever  been  reported.  The  use  of  a * S-j-Sg  is  not  unique 
here,  but  has  been  used  occasionally  especially  In  regard  to  the  effect 
of  exchange  on  spectral  llnewldths  and  lineshapes  In  concentrated  magnetic 
substances.^'14 

C.  Symmetry  Considerations 

Equations  (11)  and  (14)  are  completely  general  and  appropriate  for 
any  exchange-coupled  pair  of  metal  ions  Irrespective  of  the  relative 
orientations  of  the  single  Ion  tensors,  or  the  characteristic  single  ion 
spins  and  S2-  In  the  case  S:  * S2,  It  is  possible  that  some  crystallo- 
graphic or  approximate  symmetry  element  relates  the  two  ions  of  the  pair. 
Such  symmetry  will  In  some  way  constrain  the  form  of  both  (11)  and  (14) 
and  thus,  simplify  the  problem.  The  subject  of  symmetry  related  tensors 
has  been  considered  in  some  detail  elsewhere,15  so  will  give  only  a brief 
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description  of  how  this  bears  on  the  pair  problem.  For  what  will  follow 
we  shall  assume  that  the  single  Ion  tensors  are  real  and  symmetric;  the 
validity  of  this  assumption  Is  discussed  In  Abragam  and  Bleaney%  excellent 
monograph[16].  The  following  two  examples  should  serve  to  indicate  how  one 
goes  about  using  the  symmetry  of  the  pair  to  advantage.  Let  us  represent 
some  symmetry  related  pair  of  tensors  for  the  system  (for  the  present 
case,  the  pairs  G-|,G2,  A^,A2  and  ,D2  are  of  Interest)  as  Bj,B2.  We 
will  write  these  as  3x3  matrices  of  the  general  form 


*i 


Furthermore,  let  the  principal  values  of  B^  be  8^xx.  8izz»  and  call 

the  coordinate  system  that  makes  B^  diagonal  (the  major  axis  system  of 
B^ ) the  x.j , yi , system,  1 = 1,2.  Since  Bj  is  related  to  B2  by 
some  symmetry  element,  the  principal  values  of  both  will  be  identical, 
(l.e.,  = 8 g j * ( i , j = x,y,z)),  but  it  will  not  necessarily  be  true  that 

§ls§2  when  tw0  are  re^erred  t0  som®  common  coordinate  system.  Instead, 
they  will  be  related  by  the  transformation  equation: 

R^B-jR  = B2,  (15) 

where  R is  the  matrix  that  represents  a symmetry  element  relating  the 
pair.  We  will  now  consider  two  examples  In  order  to  Illustrate  how  (15) 
restricts  the  form  of  (11)  and  (14). 

C^O)  Symmetry 


ixx 

Bixy 

B1xz\ 

, i s 1,2. 

ixy 

Biyy 

B. 

,yz  1 

ixz 

V 

w 
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From  (15)  we  have: 


b2. 


where  i Is  the  inversion  matrix.  Without  loss  of  generality  we  may 

write  this  in  the  following  explicit  form: 

B,  B,  B, 
lxx  lxy  lxz 

Blxy  Blyy  Blyz 

Blxz  Blyz  Bizz 

Upon  performing  the  matrix  multiplication  we  have: 

Blxx  Bl*y  Blxz\ 

Blxy  8lyy  Blzz 


B2xx 

B2xy 

B2xz\ 

£ 

CM 

CO 

CO 

Sf 

B2yz  J 

(16) 

B2xz 

N 

& 

CO 

B2zz  / ‘ 

coordinate 

system,  B^.j 

' B21j 

Blxz  Blyz  Blzz/ 


(i , j , = xyz)  and  this  restricts  (11)  and  (14)  as  follows: 

(a)  The  difference  tensors  G6  and  D6  are  zero. 

(b)  The  sum  tensors  G and  D are  simply  twice  the  single  ion  tensors. 

(c)  From  (16)  it  follows  that  the  x^,y^  ,z^  system  is  identical  in  orien- 
tation to  the  x2,y2,z2  system  and  also  to  the  major  axis  system  of 

sum  tensors. 

Thus,  once  the  tensors  of  the  pair  are  determined  experimentally 
one  has  also  determined  the  single  ion  tensors  as  well;  we  will  see 
later  that  this  Is  not  always  possible.  Also,  Is  not  possible  in  this 
case  to  tell  how  the  major  axes  of  the  tensors  are  related  to  the  geometry 
of  the  pair,  nor  can  we  restrict  the  major  axes  of  to  be  the  same 
as  those  for  the  other  tensors  In  Eqs.  (11)  and  (14). 

D2(222)  Symnetry 
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For  this  case  we  have  three  mutually  perpendicular  two-folds,  two 
of  which  relate  the  two  Ions  of  the  pair,  the  third  being  collnear  with 
the  line  joining  the  two  Ions  of  the  pair.  Let  these  symmetry  axes  define 
the  x,  y,  z system  and  we  will  arbitrarily  choose  to  call  the  direction 
of  the  line  joining  the  two  ions  z.  Any  other  choice  for  this  direction 
amounts  to  a simple  relabeling  of  axes.  Eq.  (15)  now  gives  us  the  fol- 
lowing two  expressions: 


-2x  Bl-2x  = -2’  - 2y  -1  -2y  = -2* 

where  -21  represents  a rotation  of  180°  around  two-fold  parallel  to  the 
i * x,y  direction.  Performing  the  matrix  multiplication  yields: 

B2xx  B2xy  B2xz 
B2xy  B2yy  B2yz 
B2xz  B2yz  B2zz 


Blxx  ' 

Blxy  " 

BUz 

~Blxy 

Biyy 

Blyz 

-Blxz 

DO 

fsi 

Blzz 

lxx 

Blxy 

Blxz 

lxy 

Biyy 

Blyz 

lxz 

Biyz 

Blzz 

B2xx  B2xy  B2xz 
B2xy  B2yy  B2yz 

B2xz  B2yz  B2zz , 


In  order  for  both  of  the  above  to  be  true  the  following  must  also  be  true: 
Blxz  * 32xz*  °*  Biyz  * B2yz  s With  this,  the  content  of  the  above 
two  matrix  equations  can  be  expressed  as  follows, 

\ 

(17) 


Blxx  ■ Blxy 

O 

<5 

CO 

X 

X 

CVJ 

GO 

Blxy  Blyy  0 

x 

B2xy  B2yy  0 ' 

0 0 Blzz< 

1°  0 B2zz  / 
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and  (17)  Implies  then  that  z^»z2  and  z are  a1*  collnear,  that  Is,  the 
two  single  Ion  tensor  share  a common  major  axis.  Eq.  (17)  restricts  the 
sum  and  difference  tensors  to  the  following  explicit  form: 


0 0 


s* 


-2Blxy  0 \ 


2s!yy  0 


-2B.  0 

lxy 


2B 


lzz 


J 


(18) 


0 

0 0 0 

We  have  thus  obtained  a great  deal  of  information  about  how  the  major 
axes  of  the  tensors  are  related  to  the  geometry  of  the  pair.  It  Is  also 
possible  to  obtain  some  further  Information  about  the  single  ion  tensor 
themselves.  Since  both  B1  and  B2  share  z as  a major  axis,  it  follows 
that  a rotation  around  z will  be  sufficient  to  make  x^,  colinear  with 
x and  y-j  colinear  with  y.  If  we  call  the  angle  needed  to  do  this  e,  the 
corresponding  angle  for  x2  and  y2  is  «e.  In  terms  of  its  principal 
values  the  form  of  Bj  in  the  xyz  system  becomes: 


[Sfxxcos2et6lyysi"2e]  [8,xx-6lyyJs,necos9 
[Blxx-elyy]sinec°s0  [Blyycos2e+ejxxsin2e] 


The  corresponding  expression  for  B2  can  be  found  by  replacing  e by  -e 
and  1 by  2 everywhere  in  (19).  Equation  (19)  then  gives  explicit  expres- 
sions for  the  B^  and  B2i^  of  (17)  and  (18)  in  terms  of  their  principal 
values  and  the  angle  e.  Thus,  if  one  were  to  determine  experimentally 
all  the  sum  and  difference  tensors  contained  in  (11)  and  (14),  one  could 
also  find  the  principal  values  of  the  single  ion  tensors  as  well  as  the 
angle  e.  This  Is  not  always  possible,  and  to  see  this,  consider  the 
common  experimental  case  where  J is  much  greater  than  the  other  terms  In 
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the  pair  Hamiltonian.  We  thus  ignore  the  Inter-multlplet  matrix  elements 
and  use  only  Eq.  (11)  to  describe  the  pair.  Examination  of  Table  IV 
reveals  that  for  S-^Sg  the  coefficients  of  G^and  Dfi  which  are  C and  C_ 
respectively,  are  both  zero  for  all  multlplets.  It  then  follows  that 
all  the  observable  information  about  the  single  ion  G and  D tensors  Is 
contained  in  G and  0 which  from  (17),  (18)  and  (19)  have  the  following 

"CJ  “0 

form: 


B 

-o 


[B,xXc«2e*Biy/in2e] 

0 

0 


0 

[B)XJ(S1n2e*6ly)1cos2e] 

0 


Uo) 


From  (20),  it  is  clear  that  even  after  the  sum  tensors  are  found  there 
will  be  only  three  equations  for  the  four  unknowns  to  be  determined.  Thus, 
unlike  the  situation  in  C^(l)  symmetry,  determination  of  all  the  "observ- 
able" information  in  the  pair  spectra  will  not  in  general  be  sufficient 
to  completely  describe  the  single  ions  when  the  pair  symmetry  is  Ofltt). 
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Table  I 

* 

Some  Irreducible  Tensor  Operators 

Ti0(A)  - Tr(U)  = -1//2a+,  T,  _i  - 1//2a. 

T22(a)  = 1/2a l 
T21(a)  = -1/2[a2a++a+azJ 
T20(a)  = 1//6[3Az  - A-A] 

T2,-! (a)  =1/2[azA_  -A.Az] 

t2,-2(a)  = V2a2 

T22(ia)  = 1/2[z+a+] 

T2](ea)  = -1/2[zza+  + z+Az] 

T20(ia)  = 1//6[3ezaz  - z-a] 

T2>-1(za)  = 1/2[zza_  +i_az] 

T2>2(ZA)  = 1/2[Z_A_] 

T22(S1S2)  = 1/2[S1+S2+] 

T21(SiS2)  = -1/2[S1zS2++S1+S2z] 

T20<S1S2)  = 1//6C3S1zS2z-S1 *S2J 
fy(SiS2)  - 1/2 

T2-2(S1S2}  = 1/2tSl-S2-;1 

ft 

*To  find  Tjm(x)  with  x * z,  S-j,  $2  replace  a above  by  desired  operator 
variable. 


Table  III. 


Reduced  Matrix  Elements  of  T-jm(z ) and  T-jm(S^ ) 


<s,|!T1(z)||s>  = 5ss,^TsviT 


<S]||T1(S1)||S>  = 6s.s^i(Si+l) 


<S'|  |T2(z)j|S>  = 6s,s/(2S+3;(2S-1)S(S+1)/"6_ 


<Si  | |T2(S. ) | |Si>  = 6sls/(2Si+3)(2Srl)S.(S.+l)/6 


I 


(Reduced  Matrix  Elements  Continued) 


